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We discuss some important issues concerning multiphcities in quark and gluon jets in e"^e~ 
annihilation. In QCD the properties of a jet in general depends on two scales, the energy 
and virtuality of the jet. Frequently theoretical predictions apply to a situation where these 
scales coincide, while for experimental data they are often different. Thus an analysis to 
extract e.g. the asymptotic multiplicity ratio Cf/Ca between quark and gluon jets, needs a 
carefully specified jet definition, together with a calculation of nonleading corrections to the 
multiplicity evolution. 

We propose methods to systematically study the separate dependence upon the two scales 
in experimental data and compare the results with theory. We present jet finding algo- 
rithms which correspond well to the theoretically considered jets. We also show that recoil 
effects add corrections to the modified leading log approximation which are quantitatively 
important, though formally suppressed at high energies. 
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1 Introduction 



QCD predictions on the scale dependence of multiplicities in high energy qq systems are 
experimentally well confirmed [jl|. A similar multiplicity behaviour is to be expected from 
high energy gluon jets, the major difference being the different colour charges of gluons and 
quarks. In this paper we want to address a set of problems which have to be carefully treated 
for a quantitative analysis of jet properties: 

• At experimentally accessible energies subleading effects are quantitatively important. 

• The jet properties depend in general on two scales, the energy of the jet and its 
virtuality, specified by the largest possible transverse momentum of one of its subjets. 
Theoretical calculations frequently refer to a situation where the two scales coincide, 
while in experimental analyses the two scales are often different. 

• In contrast to the topology of an event, which is mostly determined by a few energetic 
particles, the multiplicities in jets are sensitive to how the softer particles are associated 
to different jets. 

As will be discussed in this paper, the multiplicity of hadrons depends not only on the 
properties of the perturbative parton cascade, but also on the soft hadronization process. 
The assumption of Local Parton Hadron Duality (LPHD) implies a direct relation between 
the number of hadrons and the number of partons, provided a proper, locally invariant, cut- 
off is imposed on the parton cascade. The asymptotic behaviour of the multiplicity in 
jets was calculated in the leading log approximation in 0]. In the modified leading log 
approximation (MLLA), subleading terms in the evolution equations of relative magnitude 
1 / Vlns are included ^ . With arguments based on "preconfinement" ^ and LPHD we 
demonstrate in this paper that an important correction factor is expected due to recoil 
effects. Although formally suppressed by a factor 1/lns, this effect is quantitatively large 
and has an essential impact on the ratio between quark and gluon jets at accessible energies. 

Many jet finding algorithms have been presented for the study of e~^e~ annihilation events 0]. 
Several of these have been successfully used in comparisons between data and theory for 
properties like the distribution in the number of jets, and how this varies with the resolution 
scale. We want to stress that our problem is a different one, as described in the third point 
above, and it will be important to specifically consider the treatment of soft particles in the 
analysis. 

The angular ordering effect in QCD implies that soft particles at large angles are emitted 
coherently from harder particles which they cannot resolve. Strictly speaking, these soft 
particles do not belong to any specific jet, and the colour factor for the emission is determined 
by the colour state of the combined unresolved partons. In ref a method (the "Cambridge" 
algorithm) is proposed to associate the above-mentioned soft particles with the quark (or 
antiquark) jet, leaving to the gluon jet only those particles directly associated with the 
emitted gluon and the gluon colour charge. In this paper we study this question further and 
propose some modified cluster algorithms. 
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If we study two-jet events obtained in e~^e~ annihilation using a jet finding algoritlim witli 
a distance measure of k±-tjpe, the jet properties (e.g. the hadron multiphcity) depend not 
only on the jet energy but also on the k±-cut used. For minimum bias events, one hemisphere 
corresponds to a quark jet where these two scales are the same (an "unbiased quark jet"). 
Similarly the multiplicity of an unbiased gluon jet corresponds to one half of an imagined 
gg system stemming from a point source. In section ^ we will discuss how this quantity is 
related to a gluon jet in a qqg event. We will also present methods by which both of the jet 
scales can be systematically examined, as well as methods designed to define unbiased jets, 
where the two scales coincide. 



Some experimental results of relevance for this paper have already been presented. In \\w\ , |1T| , 
the hemisphere opposite to two quasi-collinear heavy quark jets in e~^e~ events is analyzed. 



This corresponds well to an unbiased gluon jet, with Eg ^ a/s/2. In the energy scale 
dependence of quark and gluon jets is studied. There a fixed resolution scale is used in the 
cluster algorithm, which implies that the virtuality scale is held constant. In this paper, 
we discuss how to systematically study both scales in a jet. We also present methods to 
construct unbiased quark and gluon jets in a general three-jet topology, which enables a 
study of the scale evolution of unbiased jets. 

In our analysis we will for convenience work in the Colour Dipole Model (CDM) |13|, which 
provides a geometric picture which is easily interpreted. The outline of this paper is as fol- 
lows: In section we discuss preconfinement and LPHD. In section]^ we discuss the CDM and 
the multiplicity distributions, including corrections relevant to the MLLA approximation, 
as presented in In section ^ we discuss recoil corrections to the multiplicity evolution. 
In section ^ and § we discuss the two different scale dependences of multiplicities in jets 
and present cluster algorithms designed to extract these from data. In section |^ we present 
results obtained by MC simulations. The results are discussed in section 



2 Preconfinement and Local Parton Hadron Duality 

It is essential to realize that the hadronic multiplicity cannot be determined from perturba- 
tive QCD alone. From perturbative QCD it is possible to calculate the parton multiplicity 
within a given phase space region. As the parton multiplicity diverges for coUinear and soft 
emissions, some kind of cut-off is needed in order to find a quantity which is correlated to 
the hadronic multiplicity. The nature of this cut-off depends on the properties of the soft 
hadronization mechanism, and is therefore not calculable within perturbative theory. 

The parton cascades are dominated by planar diagrams |15[ (to leading log level in an axial 
gauge). This implies that every colour charge after a cascade has a well defined partner anti- 
charge. As shown by Amati and Veneziano 0, the final parton state can be subdivided in 
colour singlet clusters, whose masses stay limited also when the total energy becomes very 



large. In |jT6[, Marchesini, Trentadue and Veneziano showed that these "preconfinement" 
clusters are limited not only in momentum space but also in real space-time. At least in the 
large Nc limit it is then natural to assume that the clusters hadronize independently from 
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Figure 1: a) After a g^gg emission with a small angle 9, there is a screening of the new 
colour charges (ii in the picture). This implies that they give no net contribution to the 
emissions of soft gluons at polar angles larger that 9 (the angular ordering effect). Their 
contribution to the emission at smaller angles corresponds to normal dipole emission in the 
rr back-to-back frame, b) After a gluon cascade, further emissions of soft gluons corresponds 
to a set of independently emitting dipoles. 



the rest of the system. This means that the nonperturbative confinement mechanism is local 
in the sense that it combines partons which are directly colour connected and that it acts 
only locally also in momentum space. 

This idea is further developed in the notion of Local Parton Hadron Duality (LPHD), pro- 
posed by the Leningrad group |^. Here a direct relation is assumed between partons and 
hadrons, which is local in phase space, which means that the dominant features of hadron 
distributions can be obtained from the parton cascades with an appropriate cut-off which is 
locally invariant. 

A very essential feature of the QCD cascade is "soft gluon interference" and "angular order- 
ing" P]. Study e.g. a bg gluon emitting a rg gluon thereby changing its own colour to bf in a 
Lorentz frame where 9 is small (see Fig |l|a). In this case the red and antired charges screen 
each other in such a way that they do not give any emissions in directions with polar angles 
larger than 9, the angle of the first emission. For these larger angles the two gluons emit 
softer gluons coherently as a single bg gluon, while for smaller angles the two gluons emit es- 
sentially independently. (For emission angles close to 9 there is some azimuthal asymmetry.) 
This implies that in a frame where the red and antired charges move back-to-back, their 
contribution to the emission of softer gluons is just the normal emission from a separating 
charge- ant icharge pair. 

The soft gluon interference (angular ordering) and a local cut-off are incorporated in the 
Marchesini- Webber formalism for the parton cascade [|17| . This formalism is implemented in 



the HERWIG Monte Carlo |Tg, |19| . In it is also shown that a locally invariant cut-off in 
the parton virtuality is essentially equivalent to a cut-off in transverse momentum relative to 
the emitting parent parton, if this is measured in a Lorentz frame where the parton energies 
are large and the angles small. The HERWIG MC also contains a cluster fragmentation 
model with local properties. In this model the gluons are at the cut-off level split into qq 
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pairs, which are combined to colourless clusters, which finally decay into hadrons. 



A cut-ofi^ in k±_ measured in the fast-moving (small angle) frame is approximately equivalent 
to a cut-off with the same k±cut in the rest frame of the emitting charge-anticharge pair. In 
this "dipole rest frame" the emission is proportional to ^^dy, which for a given k± means 
a smooth rapidity distribution within the kinematically allowed region 



\y\ < 



where s is the squared mass of the dipole. This phase space corresponds in the fastmoving 
frame just to the region allowed by angular ordering. We see e.g. that for given k±, the 
rapidity range in the two jets is given by (for small angles 6) 



In(2j9i/A;x) -ln(cot^/2) 



In(2p2/A;±) -ln(cot^/2) 



ln(2pip2(l - cos^)/A;i) = Hs/kl), 



(2) 



i.e. exactly the same result as in the dipole rest frame. Thus in a gluonic cascade as in 
Fig 0b the emission of softer gluons corresponds to a set of independent "dipole emissions" . 
A locally invariant cut-off is obtained by a k±c\it in the individual dipole rest frames. This 
is the basis of the Lund Colour Dipole Model |]13|, which is discussed in more detail in the 
following section. 

The local feature of the hadronization mechanism is also inherent in the Lund string fragmen- 



tation model Here the hadrons which originate from the colour field stretched between 
a colour charge and its associated anticharge is independent of what happens further away in 
the system, with a correlation length corresponding to a few hadron masses ||2l|. In ref ^2 



a measure, called A, is proposed, which in the string fragmentation model is strongly corre- 
lated to the hadron multiplicity. With a local cut-off for the cascade, the A-measure is also 
strongly correlated to the parton multiplicity |2^. We note that in both cluster and string 
fragmentation there is a connection between a colour charge and its associated anticharge. 
Thus, although the parton and hadron distributions are strongly correlated, they are not 
exactly identical. 

The colour coherence and the local properties of the hadronization process are fundamental 
features of the models implemented in the MC simulation programs HERWIG 

m 



191, ARI- 



ADNE and JETSET/PYTHIA ||2^. The great phenomenological success for these 
programs in describing experimental data, in particular for e"'"e~-annihilation, is a strong 
support for the local features of the hadronization mechanisms expressed in preconfinement 
and LPHD. We note also that the independent jet fragmentation model, which does not 
have this feature, has not been able to describe the data in a satisfactory way. 

In spite of the phenomenological success mentioned above, there are still fundamental open 
questions concerning the hadronization mechanism. In the large N^. limit there is a unique 
way to connect the partons as in Fig Qd. This is, however, not the case when = 3. If two 
gluons have identical colours the confinement mechanism may connect partons which are 
not directly connected in the cascade generated by the simulation program. These "colour 
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reconnection" effects are suppressed by l/N^, and some possible consequences are discussed 



in p6|, ^ No effects have, however, been observed so far in experimental data [|lT 
Closely related to this problem is the possibility of colour reconnection between partons 
from the decay of different W:s in a W~^W~ pair at LEP2 p8|, |2^. This is of special interest 
as it might affect the W mass determination, but also here no statistically significant effects 
have yet been found, e.g. in form of different decay multiplicity or modified Bose-Einstein 
correlations 



In perturbation theory it is possible to calculate the parton multiplicity within a given region 
of phase space. The conclusion of this section is that this is not the whole story if we want to 
calculate the multiplicity of hadrons. The hadronization effects of a parton containing e.g. a 
red colour charge depends upon where in phase space its partner antired charge is located, 
and this dependence is determined by the soft hadronization mechanism. An effective cut-off 
depends on the local properties of a jet and cannot be determined as a fixed phase space 



region valid for the whole jet. In |^T| Gaffney and Mueller calculate the ratio of the parton 
multiplicity in quark and gluon jets within a fixed narrow cone, including correction terms 
of order as(Q^) (or (^(l/lnQ^)). In section | we show that, assuming a local hadronization 
mechanism based on preconfinement and LPHD, we expect a further correction term which, 
although suppressed by a factor 1/lnQ^, is numerically important. 



3 The Colour Dipole Model 
3.1 The Dipole Cascade 

A high energy qq system radiates gluons with the distribution 

_|_ 

dn = CF-^dxidx3- HtT"^ V' (3) 

2vr (l-xi)(l-X3) 



where X\ and X3 are the scaled quark and antiquark momenta, and Cp = \Nf.{\ — 



With the definitions 
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^1 = s(l - xi)(l - X3), y = hn Q _ 



. = ln(fcl/A^), !^Mm = ?^^ (4) 

271 K 



this distribution can be written 



dn = —dndy— ^ —dndy, (5) 

K 2 k 



where the last approximation holds for soft gluons. 
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Figure 2: a) The phase space for a gluon emitted from a qq dipole is a triangular region in 
the (y ,n)-plane = In k^/A"^, L = In sjls?). b) After one emission at {yi,Ki), the phase 
space for a second (softer) gluon is represented by this folded surface, c) Each emitted gluon 
increases the phase space for softer gluons. The total gluonic phase space corresponds to this 
multifaceted surface. The hadron multiplicity measure, X{L), is given by the length of the 
baseline. 




Figure 3: a) An original two-gluon system corresponds to two colour dipoles. b) After one 
emission, the emission of softer gluons (with lower kj_) corresponds to three dipoles. 



The kinematical constraint kj_ < ^/s/{2coshy) ^ ^/s ex-p(—\y\) implies that k + 2\y\ < L, 
where L is given by 

L = ln(s/A2). (6) 

Thus the allowed phase space for gluon emission is approximately a triangular region in the 
K, y-plane, cf. Fig|^. After the emission of a gluon at ki, yi, the distribution for emissions 
of softer gluons corresponds to two independently emitting dipoles, one between the quark 
and the gluon, the other between the gluon and the antiquark. The available rapidity range 
for a gluon at K2 < t^i is then ln(sqg//f;5_2) + ^^{^gq/k±2) = L + ki — 2^2- Thus the phase 
space for further emissions can be represented by a folded surface as in Fig ^d. This can 
be generalized for several emissions and a multi-gluon event corresponds to a picture with 
many folds and sub-folds as in Fig |^. 

In this language an initial gg system corresponds to two dipoles. If the gluons are e.g. rg 
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Figure 4: After one gluon emission, the different regions in the folded n-y phase space 
approximately corresponds to the angular regions shown in the right figure. The angular 
directions in which the gluon and (anti) quark emits coherently approximately corresponds to 
region B, and the emission density there is proportional to Cp. Corrections to these region 
identifications are discussed in section \6. 4 



and gr, we have a rr and a gg dipole as illustrated in Fig If one of these dipoles emits 
a gluon, the emission of softer gluons corresponds to three dipoles, as illustrated in Fig |^, 
and after n emissions we get a closed chain of n + 2 dipoles. The phase space of an original 
gg system corresponds to two triangular regions as in Fig glued together along the outer 
diagonal lines. 

The emission from a dipole stretched between two gluons is however not exactly the same 
as from a dipole stretched between a quark and an antiquark. For a gg dipole the emission 
can in analogy to Eq (j^) be described by the distribution [|1^ 

dn = — -—dxidx-i- -. (7) 

2 2tt \i-xi){l-x;) ^ ' 



For soft and coUinear emissions this goes over into the standard g— >gg splitting function . 
For soft emissions we have + x| ~ 2, which implies 

dn ~ —dndy. (8) 

This result agrees with the soft emission from a qq dipole in Eq (j^) apart from the colour 
suppressed difference between aq and ao, i.e. between Cp and Nc/2. 

After the first emission in a qq dipole, the different phase space regions for further emissions 
can be associated to different angular regions. Region A of Fig ^ roughly corresponds to 
emissions with negative rapidity in the overall CMS-frame, and region B to particles with 
positive rapidity and a larger angle to the qq direction than the first gluon. Emissions from 
region E have larger rapidity (smaller angle) than the first gluon. This is also the case for 
region C + D, with the rapidity measured in the gluon direction. The first emitted gluon 
and the (anti) quark will radiate coherently with the colour charge of the parent (anti) quark 
in region B. This argument can be generalized to a situation with several gluon emissions. 
For a cascade strongly ordered in fc^, this implies that the colour factor is Cp in the original 



7 




Figure 5: The distribution P can he subdivided into distributions for different rapidity inter- 
vals A. 

qq phase space triangle and Nc/2 on all extra folds. The identification of regions presented 
in Fig ^ is however only approximately true and we will in section study the corrections 
and their consequences in more detail. 



3.2 Multiplicity Distributions 



Assuming LPHD, the hadron multiplicity is closely related with the parton multiplicity 
n^. We will here briefly describe how the parton distribution P{n = n^, L = ln(s)) is derived 



in the dipole formulation |]T3| , in order to more easily discuss the effects of recoil corrections 
in the next section. To find the parton distribution P{n, L), we first look at the distribution 
P/\{n) in a small rapidity interval A, c.f. Fig |^. The Laplace transform 



^a(7) =^e^V>{-in)P^{n) 



(9) 



has the property 'Pai+A2 = ^Ai'^A2- Thus 

lnP(7,L)=^lnPA,(7)- 



(10) 



P/\{n) depends both on the width and the height of the interval A. We denote the phase 
space height at rapidity y by l{\y\) and define 



i?(7,/) 



lim 

A^O 



lnPA(7) 
A 



Eq (1^) can then be written 
lnP(7,L): 



dyR{iA\y\)) 



^max 



In [22 1, it is shown that R also is related to V by 



^^-i?«(7,L) = ^[pfe)(7,L) 



2 / dlR{-f,l) 



q, g, «g 



dy 



dl 



"0- 



111 



(12) 



(13) 
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Figure 6: The hyperbolic shape of the true phase space limits and the inequality + < 
+ X3 < 2 are both neglected in the leading order result. This can be corrected for by 
cutting off a strip at the triangle edges, thus reducing the available phase space. The different 
heights of the strips reflects the difference between + x^ and x\ + x1, appearing in the 
emission density for a gg- and qq dipole, respectively. The different magnitudes of phase 
space reductions implies n('i) (L) ~ nfe) (L + Cg — Cq) . 

The boundary condition at some cut-off scale for the cascade 

P^i\n,ti,)=5r,o R^'\i,K,) = 0, (14) 



then imphes 



(15) 



Combining Eq ([T^) and Eq (|l^) gives the following differential equation for V, valid in the 
LLA: 



^^2-lnP«(7,L) = -^[pfe)(^^2.)-l 



3.3 MLLA Corrections 



Since the approximate triangular (k, y) region is somewhat larger than the true hyperbolic 
shape of the phase space, and the inequality xl + x^ < xf + x"^ < 2 was neglected, the 
emission density is overestimated in LLA (especially in a gg-dipole). In the Modified Leading 
Logarithmic approximation (MLLA) M, corrections of relative order 1/vT are included. 



In IjT^ it is shown, that evolution equations correct to relative order 1/vL are obtained if 
we maintain the approximation + Xg ~ + x| ~ 2 and instead reduce the available phase 
space by cutting out a strip at the edges, as illustrated in Fig ||. The height of the strip is 

Cq = 3/2, Cg = 11/6, (17) 

for a qq- and gg-dipole, respectively. The constant phase space reduction modifies |d?//d/| 
to 

2|d2//d/|« = e(L-Q-0, (18) 
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which imphes 



^lnP»(7,L) = a,i?«(7,L-Q), (19) 
'^'rlnP«(7,L) = -^[pfe)(7,L-c,)-ll. (20) 



dL2 ^" ^ L-c- 



For the gluon case, the equations are modified by the possibihty of g-^qq sphttings |T^. We 
do not reproduce the algebraic details here, but refer to refs [H, Extracting moments in 
7 of Eq (|20| ) leads to differential equations for n^^ and n^^). Their asymptotic behaviours 
are 

ttH ^ {aja^yd^ ~ LP exp(2y«^), (21) 

where 

P = ^-(2iVf/iVe' + ll)^ (22) 



in MLLA. Thus the introduction of 1/ vL-suppressed terms in the evolution equations 
changes the asymptotic behaviour of n('). 

The contribution to p from the g— *>qq process is the term 2Ni/N^ in the parenthesis, which 
is clearly a small contribution compared to the other term 11. The contribution from this 
process to the ratio N^/Ng is also very small. Using numerical calculations we have found 
that it only modifies the ratio by less than 2% for all energies above 4 GeV. For this reason 
we will neglect the splitting process in the analytical calculations presented below. However, 
when we compare our analytic results with MC simulations, the latter will include also the 
g^qq process. 

Neglecting the g^qq process for these reasons, we find from Eq ( PP]) for the ratio N^/Ng the 

MLLA result 

— n(q) (L) ^ nfe) (L + Cg- Cq) , (23) 

which essentially is a correction approaching the LLA result Eq (|2ll) for large L. 

We also want to point out that the distribution in the multiplicity measure A, mentioned in 



section |3.2| , satisfy exactly the same evolution equations as the distribution in n. Therefore 
the asymptotic increase is also the same. The boundary conditions at threshold are however 
different, which implies some deviations at lower energies. 



3.4 Virtuality Scale Dependence 



If we consider a two-jet event sample, picked from qq events using some jet resolution scale 
k±r, the mean multiplicity of the sample will be related to phase space area below = 
2ln{k±r/A). As shown in Fig|^, this implies 



n 



(^)(L, ft < Kr) 



n 



Mi 



+ (L-ft, -Cq) ^n(q)(/ + Cq) 



(24) 
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Figure 7: a) A two-scale dependent multiplicity can be studied in two-jet qq events, using 
a resolution 7^ L. The multiplicity is related to the phase space below the cut-off line 
at Kr- The parton multiplicity is given by ?7,Wq,+^ = ?7,(')(fi;r + Cj), and = {L — Kr — 
Cj)^nW(Kr + Q). b) If the hardest gluon jet is found at Kr, the phase space which determines 
the multiplicity is given by the same region a + /3 + 7, plus a fold. The multiplicity on this 
fold is given by nfe)(Kr). 

Here the first term on the right hand side corresponds to the region a + (3 m Fig while 
the second term corresponds to the region 7. Thus the mean hadron muhiphcity, A^'^, in qq 
events where no jet is resolved above is given by 

iVq\(L, K<Hr) = Nl^^{Kr + Cq) + (L - ft, - Cq)^iVq\(ft, + Cq). (25) 

Events where the hardest gluon jet is found at corresponds to the phase space below the 
horizontal line in Fig|^. The hard gluon jet gives the contribution }^N^^[Kr) and therefore 

N^^{L, K = Kr)= N^^{L, K<Kr) + ^N^^i^r) ■ (26) 

We note that while the Kr dependence is fairly complicated, the L dependence is simply 
linear. 



4 Recoil Effects and Boundary Conditions 

The ratio of the parton multiplicities in quark and gluon jets within a narrow cone was 
calculated including corrections to order a^, i.e. of order 1/L, in ref (c.f. also ref |l32| ). 
The result is that the 0{1/L) corrections are very small, only a few percent for energies 
above 30 Gev. 

In a parton cascade, A cannot be associated to a specific renormalization scheme, but must 
be treated as a free parameter. Changing A introduces terms of relative suppression 1/L 
in the evolution equations Eq (pOf). Within the language of parton cascades, it is therefore 
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non-trivial to systematically find a complete set of 1/L suppressed correction terms to many 
observables. We note however that the 1/ L terms induced by changing A are similar for the 
multiplicity in quark and gluon jets. Therefore these terms cancel to a large extent for the 
ratio Nq/Ng. This is in agreement with the result in ref |3ll], that the correction term of 
order to the ratio Nq/Ng is small and scheme independent. 

As discussed in detail in section ^, the observable hadrons are not directly determined by the 
partons within a fixed phase space region. The hadrons originating from a particular parton 
colour charge depend also on the location in phase space of its corresponding anticharge. 
This implies that the necessary cut-off in the parton cascade should depend upon the location 
in phase space of an associated colour anticharge. Such a dependence appears naturally in 
the dipole cascade formalism with a fc^-cut in the dipole restframe. Therefore the result 
can be sensitive to recoil effects experienced by the emitting parent partons. Assuming the 
local properties of the hadronization mechanism discussed in section 2, we will in the next 
subsection derive a correction term to the ratio N^/Ng. Although formally of order 1/L, 
this recoil effect is numerically large, and important for comparisons between theory and 
experimental results. 



In section we compare our results with data from OPAL [TT] for jets with energy Mz/2. 
In section ^ we also discuss how to analyze LEP data at the Z pole in order to obtain the 
ratio Nq^/Ng for a range of jet energies. 



4.1 Recoil Effects 



As discussed in section O we will neglect the effect caused by the g— *>qq process. Numerical 



calculations have shown that this process, which is both colour suppressed and kinematically 
suppressed, only has a negligible influence on the result. In principle the splitting process 
could be included at the expense of more complicated expressions, which would make the 
result less transparent. 



Let us study the emission from an original gg-system. As discussed in section |3.1| this 
emission corresponds to two dipoles (a rr and a gg dipole in Fig ^), and the phase space 
can be represented by two parallel triangular regions. If one gluon is emitted (e.g. from the 
gg dipole), the ensuing 3g state emits further soft gluons as three dipoles, as indicated in 
Fig 0b. As discussed in section 0, the emission from the separation between the b charge in 
the emitted gluon (2) and the b charge in the emitting (recoiling) gluon (1) is determined by 
the invariant mass Si2. The Landau-Pomeranchuk formation time for an emitted quantum is 
determined by the k± of the emission. Thus there is a time ordering, such that a softer gluon 
(a gluon with smaller k±) is emitted after the establishment of the configuration in Fig ^d. 
We must then expect that the weight of the rr and gg dipoles in Fig are determined by 
the corresponding dipole masses. The squared mass of the "spectator" rf dipole is reduced 
from its initial value s to sis = (1 — X2)s, which thus reduces the phase space for emissions 
from this dipole. We note that this recoil effect is very different for a qq system, where the 
weight of the corresponding dipole is negative, and suppressed by a factor Therefore 
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recoils will affect the ratio between quark and gluon jets, and we will now estimate the size 
of this effect. 



The density of gluons with squared transverse momentum k'j_ = s{l — Xi){l — x^) emitted 
by e.g. the gg dipole in Fig ^a is given by Eq (|^). This initial emission reduces the emission 
of softer gluons with k'l < k\ in the partner rr dipole, so that the rapidity range is reduced 
from Ay = Ins /k'l to Ay = Ins 13 /k'l = Ins /k'l + ln(l — X2). On average the reduction 6y 
caused by gluons within an interval dln/c^. = d/t is then 

5y = dK / dxidxs— ^ , rln(l - X2)5iK - ln[s(l - Xi)(l - X3)]) = 

K J 2(1 — — X3) 

= ^dK-J(L-K), {L = \ns). (27) 
For a strongly ordered cascade we have L ^ k, and in this limit we find 

ECr, L-K>1. (28) 




Deviations from this limiting value give corrections which are suppressed at high energies. 
We have however also checked that the value in Eq (^) is a very good approximation 
in the whole relevant part of phase space. Typical values of L — k, become smaller for 
smaller L-values, i.e. smaller dipole energies, due to the increase of as- Within the MLLA 
approximation we find that for L ~ 6 the average value of k becomes ~ 2, i.e. k± ^ eA, 
which is very close to the standard cutoff in the ARIADNE cascade MC. Thus typical values 
of L — K, are always larger than 4, and for these values / deviates from the value in Eq (|28| ) 
by less than 10%. 



Half of the effect in Eq (|27D corresponds to positive rapidities, in which case the recoil is 
taken essentially by the gluon (1) in Fig |^, while in the other half the recoil is taken by 
gluon (3). The emissions in the rr dipole in Fig ^ has however a similar effect on the gg 
system, and the net result on one of the gluon jets is thus given by Eq (p^). 

Once one gluon is emitted two smaller dipoles are formed (between the pairs 3-2 and 2-1 in 
Fig Pd) emitting softer gluons, which give similar reductions of the "spectator dipole" phase 
space. Further down the emission cascade the original gg dipole in Fig ||a corresponds to 
a chain of dipoles. Emissions in the two end links of this chain give recoil effects on the 
spectator system originating from the rr dipole in Fig Consequently we obtain a recoil 
effect described by Eq (^) for all values of k, not only those which correspond to the first 
emission. Summing up the effects from all emissions above a given value / = In/c^/A^, and 
replacing / with its limiting value Cj., gives the total phase space reduction 

Sytot{l) = / dn — Cr = aoCrln — z-^. (29) 

The upper limit is given by the kinematical limit in the MLLA approximation, L — Cg. Thus 
the effective rapidity range becomes 

L — c 

Ay^sil) = L - I - - aocM — — ^, (30) 
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where Cg is the MLLA correction and the last term is the recoil effect. We note in particular 
that the recoil term goes to zero at the kinematical limit / = L — Cg, as it should. From 
Eq m) we find 



dl 



^, (31) 



which should be inserted in into Eq (|^) (with 2y replaced by A?/efr)- Doing so, and taking 
the MLLA modification in Eq (|TB]) into account, we find instead of Eq (|TP|) 



d 

dZ 



InPfe) (7, L)=(^l- (7, L - 



Cg). (32) 



We note in particular that the result in Eq (^) is independent of L and that the correction 
in Eq (|29| ) goes to zero in the kinematical limit. These two features imply that the structure 
of the evolution equation is unchanged when going from Eq (pj]) to Eq (|32D . 



For an initial qq system the emission from a quark colour charge is lower than one half of a 
gluon charge, Cp = — For the radiation from a qqg system, this can be expressed 

as a negative contribution with relative weight from a dipole stretched between the 

quark and antiquark, corresponding to the rr dipole in Fig for the purely gluonic case. If 
we assume that we can treat this negative dipole in the same way, we would get the following 
result for a quark jet 



Alnp(.)(^,L)=(l + -^)i?('^)(7,L-c,), 



(33) 



- -4 t dxi±ii^ln(l - X) = ^ ( ^ - ^ 1 . (34) 



iV2 Jo X 

Combining Eqs (p2|), (B3) and (pIS]) now leads to 




dL aqi + aocr'/{L -Cq) dL 

« ^(l-^^)^toP'"(7.L). (36) 

where Cr now is modified by approximately 10%, to 

Cr ^ + = f^n^-l (36) 



The negative dipole in the qqg system is caused by interference effects due to identical 
colour charges. As discussed in section 2, these interference effects are also connected to the 
possibility of "colour reconnection" in the hadronization process. As these problems are still 
not solved, we regard the recoil effect for the quark jet as uncertain. Since this contribution 
is only 10% of the total effect this is however not a serious problem here. 

Taking the first moment in the Laplace transform variable 7, we find the relation Eq (P5|) 
between jfn^^^ and -^n(^) as well. Remembering that rife) refers to a gluon jet while n'^^^ 
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refers to a qq system, the relations between multiplicities in two-parton systems should be 
given by 



4.2 Boundary Conditions for Hadron Multiplicities 



The relation Eq ( pTl) has to be supplemented by appropriate boundary conditions. Extrap- 
olating Eq (|37|) down to too small values of L would imply that the hadron multiplicity in 
a qq system would be significantly larger than that in a gg system. At low values of L, 
the hadron multiplicity is largely determined by the hadronic phase space and thus by the 
total available energy. This implies that at some threshold value Lq, we expect to have the 
relation 

iVg^g(Lo) ^ iVq\(Lo) = iVo, (38) 
Here Lq, though larger than Kc-, should correspond to an energy of only a few GeV. 

The precise value of Lq strongly depends on non-perturbative QCD effects. At low energies, 
the fact that a qq string contains two quarks while a gg string does not, may influence 
the ratio N^J N^^. Thus the value of Lq is sensitive to details in the fragmentation of 
low energy qq and gg systems, while A^o to a large extent depends on how the primarily 
produced hadrons decay. In principle Lq ought to be determined by experimental data 
from charmonium and bottonium decays. In our analysis, we have instead determined Lq 
and A'o from Monte Carlo simulations of the Lund String Fragmentation model, using the 



JETSET 7.4 computer program P5|. We then get Lq ~ 5.7, corresponding to an energy 



Ea = Aexp(Lo/2) ~ 4GeV for A = 0.22GeV. Given N^^ and Lq from MC simulations, N^^ 
can be derived by numerical integration of the right hand side of Eq (0). The dependence 
on the threshold behaviour can however be avoided if one studies how the multiplicity varies 
with increasing energy. This possibility is further discussed in section 0. 



5 Scale Dependences in Jets 



In the following discussion we will use the notation of and i.e. and Nj^ denote 



multiplicities in jets, while A'gg and refer to multiplicities in two-parton systems. We 
have already introduced L and k as logarithmic energy- and virtuality scales for two-parton 
systems. The multiplicity in a jet j with energy Ej will be studied as a function of the 
logarithmic jet energy scale 

L, ^ In {{2E,r/A') . (39) 

The jet energy is multiplied by a factor 2 for two "cosmetic" reasons: Then the scales coincide 
for a gluon jet when the simple condition 

L, = K (40) 
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Figure 8: The transverse resolution scale k± between the softer jets in a three-jet event also 
specifies the maximal allowed transverse momentum of unresolved particles within the jets. 
The multiplicity of the jets thus depends both on the jet energy and on k±. 



is fulfilled. (I.e. the relevant energy scale of this gluon jet is Eg = k±/2, which is the MLLA 
result presented in |^.) Furthermore, the multiplicity in a one-scale dependent jet is simply 
given by 

iVp^(L,) = V (L,), p = q,g. (41) 



In general, the multiplicity of a jet depends on two scales, energy and virtuality. Reducing 
the maximal allowed transverse momentum within a jet reduces the multiplicity, even if the 
energy remains constant (c.f. Fig |]). 



As discussed in section the virtuality scale dependence of the multiplicity in quark jets 
can be easily studied in two-jet qq events varying the resolution scale. To compare quark and 
gluon jets we can use a three-jet event sample. The multiplicity in the jets will also in this 
sample depend both on the jet energy and the resolution scale k±r- Keeping k±r fixed, the 
energy scale dependence can be studied |T^. However, the restriction implied by a fixed kj_r 
is rarely considered in theoretical predictions. Rather, most calculations apply to the very 
forward region of jets |^, |3l|]. There coherence effects are negligible, and energy-momentum 
conservation restricts k± of emissions to such an extent that k±r introduces no new bias. 
We will here examine the possibility to study unbiased jets at moderate energies, using the 
broadest possible cones to define the contents of the jets. 

To find the jets of an event, cluster algorithms are used. In general, these contain a definition 
of a distance measure d between the jets, combine the two closest jets into one, and continues 
until all distances are above some resolution scale dcut- In CDM and other approaches [|l^, 
the transverse momentum specifies the resolution. Thus it is appropriate to use a k^-tjpe of 



distance measure, as e.g. in the Durham LUCLUS p5| or DICLUS algorithms. We 
will use an approach, where jets are combined in the order given by the chosen algorithm 
until only three remain. The topology of the jets then specifies the transverse momentum, 
k±, of the event (c.f. Fig |^). This approach enables the definition of unbiased jets, whose 
evolution with k± can be studied, while the energy scale dependence can be examined in an 
event sample with fixed k±. 

To construct unbiased jets, we note that a cone-like surface used to define the forward 
contents of a jet corresponds to a perpendicular plane in some other frame. In this frame 
the jet resembles one hemisphere of a full event, and if the energy and virtuality scales in 
this frame are equal, the jet is unbiased. We will refer to this observation as the "One-Scale 
Criterion" (OSC), since the mean multiplicity of the jet in this case indeed only depends on 
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one scale - the energy of the corresponding full event. 



6 Jet Algorithms 

Soft gluon coherence can be approximated by an angular ordering (AO) constraint |p. The 
multiplicity in a jet defined by AO will thus depend on one colour factor only. According to 
LLA, the jet will also be unbiased and depend on one single transverse energy scale, provided 
no extra cut-off in transverse momentum is imposed by a fixed jet resolution scale. We will 
present results from two algorithms based on this AO observation: The Cambridge algorithm 
presented in ^ , and a previously undiscussed "Mercedes" algorithm, where the multiplicities 
of the jets are defined in the symmetric "Mercedes" Lorentz frame of the event. The main 
reason for studying both is to look for similarities, which point at general properties in the 
AO approximation. 

Corrections to the AO algorithms can be found by using the OSC explicitly. We will here 
discuss one such algorithm, which has similarities with the Cambridge algorithm, called 
the "Cone Exclusion" algorithm, and one corrected "Boost" algorithm where the jets are 
analyzed in a frame similar to the Mercedes one. The two OSC algorithms give the same 
gluon jets, but the definition of quark jets differ. With the Cone Exclusion algorithm all 
three jets are unbiased, while the Boost algorithm constructs quark jets with well-defined 
but different energy and virtuality scales. Thus the Boost algorithm is suitable for a study 
of the separate dependence of the two scales, while the Cone Exclusion algorithm is suitable 
for a comparison of quark jets vs. gluon jets. 



6.1 Angular Ordering Algorithms 

The Cambridge Algorithm 

The Cambridge algorithm is designed to construct gluon jets uncontaminated by coherently 
emitted particles. Thus the gluon jet properties depend only on the transverse momentum 
to the nearest harder jet, and on the gluon colour factor Nc- More specifically, the Durham 
/c_L-distance 

d.,^2 LJ-^(1-cos%) (42) 

is used to resolve jets, but the particles and sub-jets are merged in inverse angular order 
(those closest in angle are combined first). Once a soft jet is resolved, it is "frozen out", i.e. 
it gets no extra multiplicity contribution. Thus the contents of the soft jet is confined to 
a cone given by the smallest angle to any harder jet. This implies that the treatment of a 
gluon jet depends on wether its energy is higher or smaller than the quark jet energy. The 
Cambridge algorithm is therefore suitable for gluon jet analyses primarily when the hardest 
gluon jet is significantly softer than the quark jets. 
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In the Cambridge algorithm, dcut not only determines when the clustering stops, but also at 
what stage different jets are frozen out. Producing a fixed number of jets is therefore not 
a completely trivial task. If we put d^nt artificially large and pursue the clustering down to 
three jets, no freezing will occur. This would then correspond to a strict angular ordered 
clustering, similar to a cone algorithm. Fixing the number of jets is therefore better achieved 
by changing dcut in every event to a value which produces three jets. This procedure faces 
two problems: There may not be any dcut giving three jets, and when there is a large range 
of dcut values giving three jets, the multiplicities in the jets could depend on the choice of 
dcut- However, both of these situations occur very rarely (at percentage level), and reliable 
conclusions may therefore be drawn from events clustered by the Cambridge algorithm, 
varying dcut to fix the number of jets to three. 

The Mercedes Algorithm 

In conventional cluster algorithms the bisectors between jets roughly separate the contents 
of them. In the specific case of a completely symmetric three-jet event, commonly referred to 
as a Mercedes event, such a jet definition will satisfy angular ordering [Q. The kinematical 
constraint implies a fixed scale, but the scale evolution can still be studied by boosting a 
general three-jet event to its Mercedes frame. Since particles are in general shuffled from one 
jet to another under a Lorentz transformation, the mass of a jet is not invariant. Instead 
the direction of a jet, corresponding to a parton in the cascade, approximately transforms 
as a light-like vector. Using a /c^-based algorithm to find the jets, the Mercedes algorithm 
constructs gluon jets similar to those in the Cambridge algorithm. 

6.2 OSC Algorithms 

The Cone Exclusion Algorithm 

The "Cone Exclusion" (CE) method combines k±- and angular distances in a way that has 
similarities with the Cambridge algorithm. After the construction of three jets using a k±- 
distance, a cone-like region is defined around the gluon jet. Only particles assigned to the 
jet that lie inside the region are then contributing to the multiplicity. Note that in spite of 
the name "Cone Exclusion", the method is using a fc^-based cluster scheme to find the jets. 
The "Cones" are used to assign soft particles to the correct jet. 

With this method it is simpler to fix the number of jets to three than in the Cambridge 
algorithm, and it provides a better treatment of hard gluon jets. It is however specifically 
designed for studies of multiplicities in jets in three-jet events, and does not share the benefits 
of the Cambridge algorithm as compared to other fc^-algorithms in other respects. E.g., the 
Cambridge algorithm is designed to avoid the formation of "junk jets" (when soft particles 
from different jets which happen to be close in phase space are combined and may be resolved 
as a jet if the resolution scale dcut is small). Since we in this analysis always pursue the 
clustering until only three jets remain, these "junk jets" will in general be absorbed into 
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Figure 9: a) // the partons of a dipole are going apart back-to-back, particles from the trian- 
gular region of height Kr in the direction of jet j all have a rapidity y > \ia(2Ej/A) — Kr/2, 
which implies cot{9j/2) > 2Ej/k±r- This holds also if the dipole is not in its CMS. b) If the 
dipole is in a general frame, the triangular phase space areas corresponds to two "egg-shaped" 
regions which after a boost along the bisector to the back-to-back frame become the cones of 
(a), i.e. which satisfy cot{6j/2) > 2Ej/kj_r, where 6j and Ej are measured in the collinear 
frame. 

"proper jets", and we are therefore not very sensitive to this problem addressed by the 
Cambridge algorithm. 

To find the proper cones to use, we will start our discussion with two-jets events obtained 
with some resolution k±r. We look at the events in a frame where the original partons are 
going out back-to-back, though not necessarily in the overall CMS. In the K,y phase space 
picture (with folds), an unbiased jet corresponds to a triangular region. When the partons 
of a dipole are moving apart back-to-back, the relation AE1E2 = s implies 

9 7? 9 /? 1 

L = H^) + \n{^) = -{L, + L,), (43) 

where Li and L2 are the logarithmic jet energy scales in this frame. A particle p stemming 
from a triangular region of height = 2ln{k±r/-^) in the direction of e.g. jet 2 has an angle 
9p2 to the jet given by (c.f. Fig ^) 
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An equivalent expression for this constraint is 

In(spi) - ln(sp2) > ln(si2) - ln(A;^J. 



(45) 



We now turn to three-jet events, where k\ = s(l — xi)(l — X3) specifies the resolution. Then 
all quantities in the constraint in Eq (^5]) are Lorentz invariant, and it can immediately be 
applied to the three jets in the CMS. For the qg dipole, the constraints are then 



for the gluon jet and 



Xg sin^(gpg/2) 
Xq sin^(6'pq/2) 

Xq sin2(6'pq/2) 



< 1 - Xq (46) 



< 1 - Xq (47) 



Xg sin-(^pg/2) 

for the quark jet. Here ^pq denote as before the angles between the hadron p and the jet 
direction q or g. The relations for the qg dipole are similar. 

For the q (or q) jet, Eq (|^) specifies an "egg-shaped" region, which after a boost along 
the bisector of the dipole to a back-to-back frame becomes a cone satisfying Eq (^ij) (c.f. 
FigP). The gluon jet is however attached to two dipoles, and CE thus specifies two different 
regions. In the directions being accepted by one dipole but not the other, most particles 
emerge from the "wrong" dipole and should not contribute to the multiplicity of the gluon 
jet. We have required particles belonging to the gluon jet to satisfy both possible restrictions 
from Eq 



To summarize, the Cone Exclusion algorithm works as follows: Three jets are constructed, 
using a k±-hased cluster algorithm. A particle assigned to the gluon jet will then contribute 
to the multiplicity only if it satisfies Eq (0). The quark jets are treated similarly, using 
Eq (^Tf) . Thus the multiplicity in the forward region of every jet - corresponding to an 
unbiased jet - is studied, while soft central particles are simply ignored. 



The Boost Algorithm 

The OSC can also be used to improve the AO-based Lorentz transformation to the Mercedes 
frame. Consider a three-jet event boosted to a Lorentz frame where the angles are 

^gg = ^qg = (48) 

and the partons carry energies E'^. Let the bisectors define the planes between different jet 
regions. For the gluon, the logarithmic back-to-back jet energy scale from both dipoles is 
then 

2E' 9' 

Lg = 21n(-^sm-), (49) 
which coincides with the virtuality scale when Lg = k. This implies 

x^^sin2(072) = (1-Xq)(l-Xq), x'i = ^. (50) 
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The search for cone-hke boundaries which define an unbiased gluon jet has thus been refor- 
mulated to a search for a specific Lorentz frame where the jet regions are easily identified. 

After a bit of algebra, one finds the general relation 

x'^sm^{e'/2) = (l-a;q)(l-Xq)-^— cos2(072). (51) 

1 — 

Thus the requirement Eq (|50|) is satisfied when 

cos^(^72) = (52) 

In the soft gluon limit, the wanted Lorentz frame coincides with the Mercedes frame, with 
6' = 120°, but larger gluon energies give larger angles 6'. We also note that the energy scale 
for the quark jet in this frame is 

<sin^(^72) = ^, (53) 

and similarly for the q jet. 

The gluon jet defined in this way is actually equivalent to the one defined by the CE algo- 
rithm. If we combine Eq ( ^O]) and (^) into x'Jx'^ = 1 — Xq and exploit the Lorentz invariance 
of the left hand side of Eq (^Bj), we note that the CE condition for the gluon jet in Eq ( P6| ) 
can be rewritten as sin^(6'pg/2) < sin^(6'pq/2), which is just the bisector condition in the 
Boost algorithm. The treatment of massive particles, whose masses have been neglected 
in the discussion, differ however in the two algorithms. Using both is a simple test of the 
sensitivity on particle masses. 

The definition of quark jets differ in the two OSC algorithms. From Eq ([53|), we see that 
the quark jets in the Boost algorithm correspond to regions A^B and E in Fig ^. Thus the 
Boost algorithm provides means to study the two-scale dependent multiplicity for different 
jet energy scales in a fixed energy experiment, while the CE algorithm is better to use for 
the study of one-scale dependent multiplicities. 

To summarize, the Boost algorithm is as follows: Find three jets using a fc^-based cluster 
algorithm. Boost the event to the frame where the jet directions, assumed to be massless, 
satisfies Eq ( ^8] ) and Eq (|^). Let the jet boundaries be given by the bisectors to the other 
jets in this new frame and re- assign particles to the jets accordingly. 



7 Results 



We have tested the analytic form for the ratio Nj^^/ N^g in Eq (^) by comparing with 
multiplicity results from MC simulations and preliminary data |T^. Using MC simulations, 
we also examine the presented jet algorithms by comparing the multiplicities obtained in 
jets with the multiplicity of complete events at corresponding energies. 
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The CDM is available as a Monte Carlo simulation program, ARIADNE [2^]. There dis- 
tribution factors such as xl + x\ are taken into account and energy conservation is obeyed. 
Thus Monte Carlo simulations can show whether corrections other than presented above are 
needed to understand the model predictions on the scale evolution. 

As discussed in section ^, a qqg system can be regarded as three dipoles, where the qq dipole 
is colour suppressed and has a negative weight. Alternatively, the system can be described 
by two dipoles, qg and gq, where the colour factor transforms from Cp in the q (q) end to 
Nc/2 in the gluon direction. However, dipoles with negative or non-uniform colour factors 
are ill suited for MC implementation. In the standard ARIADNE Monte Carlo, the solution 
is to neglect most terms of order A*""^, using the colour factor Nc/2 in all gg- and qg dipoles. 



In a modification to the MC correcting for this approximation is presented. There 
non-uniform colour factors in dipoles are implemented, in a way reflecting the discussion 
around Fig ^, where the emission density is assumed to be proportional to C-p in all of the 
original qq phase space triangle, while Nc/2 applies to all extra phase space folds after gluon 
emissions. 

At moderate energies, the mean multiplicity is mostly determined by the hardest gluon. The 
emission of this gluon in the qq dipole is correctly given by C-p in both MC versions. Thus, 
results from the two MC approaches are expected to deviate only at larger energies, showing 
a relative discrepancy of at most 

A ggg configuration is well described by three dipoles, all with positive emission density, 
determined by Nc/2. This picture is implemented in default ARIADNE, which is used to 
obtain the simulation results from gg-systems presented below. 

In the simulations we have used the default tune of ARIADNE 4.08, with A = 0.22GeV. The 
parton configurations obtained in the cascade simulation are hadronized using the JETSET 
Monte Carlo [53], which is an implementation of the Lund String Fragmentation model. 



7.1 The Multiplicity Ratio N^jN^g 

From MC results we find the threshold energy where the multiplicity in quark and gluon 
jets are the same, to be given by Lq ~ 5.7. This corresponds to a CMS energy of ~ 4GeV. 
Lq also specifies A^^q = N^^{Lo). 

In Fig |TU|a, the simulated ratio of N^^^/N^^ is compared to the prediction of Eq (^7\), using 
values of A'o, Lq and N^^{L) obtained from the MC. If the 1/L recoil correction in Eq (|37|) is 
neglected the result differs significantly, even at L ~ 20, i.e. at y/s ~ 5TeV. The prediction 
is further modified if the MLLA correction of order 1/vT in Eq (pTf ) is neglected. 

Fig |TU|b shows the different multiplicity ratios obtained in simulations of e^e~ annihilation 
and pure dd events. At energies slightly above a heavy quark thresh-hold, the mean mul- 
tiplicity gets a significant contribution from isotropically decaying heavy hadrons. This is 
seen as a peak at L ~ 6 {^/s ~ 2m£)) and a little shoulder starting at L ~ 8 {^/s ~ 21713)- 
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Figure 10: Multiplicities ratios o/qq and gg systems, based on all stable particles, qq samples 
are based on pure dd events (diamonds) and e^e~ annihilation, using ARIADNE default 
(boxes) and the modification of ^Ij] (crosses). The statistical errors for the MC results are 
within symbol sizes. We find that N^^^ = N^^ at Lq = 5.7. a) Comparison with the prediction 
of Eq ^3li), using values of Nq, Lq and N^^{L) obtained from the MC. If neglecting the recoil 
correction (dashed line) and also the scale shift Cg — Cq (dotted line), the prediction is far 
from MC. Thus recoil corrections of order 1/L are essential, even at ^/s ~ 5TeV. b) For 
e'^e" events, cc and bb threshold effects are seen as a peak at L 6 and a shoulder starting 
at L ^ 8. The relation determined by Eq j ^Slj ) fits well with Lq = 6.7. The difference between 
default and modified MC are important at L 20, but negligible at Zq energies and below, c) 
The -^N{L) ratio approaches the asymptotic value Cp/Nc faster, but subleading corrections 
cannot be neglected, d) The ratio of multiplicity derivatives -^N^L) is very similar for dd 
and e~^e~ simulations and thus independent of Lq, as expected from Eq l^3l{). 



The details of these threshold effects are beyond the scope of this paper. Instead we note 

6.7 (v/i 



that Eq (0) still fits well with Lq 



the corrections to N^^^ from the Monte Carlo modification of 127 



6.3GeV). In Fig |T0|b it is also shown that 
which more consequently 



implements the difference between Cp and Nc/2 after the first gluon emission, are important 
at very high L, but negligible at Zq energies (L = 12) and below. 



The ratio between the derivatives j^N{L) is presented in Fig p!0|c. This ratio is expected 
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to approach the asymptotic value Cp/Nc more rapidly, which is also confirmed by data [12]. 
This is also born out in our analysis, and we note a good agreement between MC and the 
analytic form in Eq (p?]). According to this relation, the ratio of multiplicity derivatives 
-^N(L) is expected to be independent of Lq. Fig |lO|d shows the similar results of dd and 
e~^e~ simulations, confirming this expectation. 

In |Tl|, the OPAL collaboration has studied the /Nj^ ratio via quark and gluon hemi- 



spheres with energies E = Mz/2 [L = 12). Their preliminary result 



qq 



(L = 12) = 1.509 ± 0.022(stat) ± 0.046(s?/st) (54) 



is in excellent agreement with our result including recoil effects in Fig p!0|b, 

g(i.l2,.0.67^^. (55) 

In |Tl|, OPAL also present preliminary data for the ratio of multiplicities in the central 
rapidity region, 

-r^{L = 12, \y\ < 2) = 1.815 ± 0.038(stat) ± 0.062(s?/st). (56) 
qq 

For a reasonably small central rapidity range Ay, the multiplicity N{L, Ay) corresponds well 
to Ayj^N{L). It is therefore interesting to compare this result with the analytical expression 
for the ratio of multiplicity derivatives, which is more independent of the boundary conditions 
at Lq. From Fig |TO|c, we find 

^^^(L = 12) = 0.547 = (57) 

with recoil effects included. The MLLA prediction without recoil effects is smaller than 1/2. 

Thus we conclude that MLLA calculations complemented with the recoil effect discussed in 
this paper are in very good agreement with the preliminary experimental data from OPAL. 
The recoil effect gives a sizeable correction which implies that the asymptotic value Cp/Nc 
is far beyond reach in accelerator experiments. 



7.2 Determination of Multiplicities in Jets 

The quantity N^^ in the previous subsection is the multiplicity in a gg event which is hard 
to realize and study directly in experiments. In section ^, different jet algorithms designed 
to investigate this quantity in normal three-jet events are presented. We have tested their 
performance by simulating events at Zq energies. The events are clustered into three jets 
and the multiplicities in the jets from different algorithms are studied as a function of k and 
y. The results for the jets are compared with the full-event results presented above. 
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Figure 11: a) Multiplicities in gluon jets as a function of k. Results for several values of 

y are plotted, or the range of results are represented with shaded areas. The AO algorithms 
(Cambridge and Mercedes) have similar behaviours, and the multiplicity approaches O.bN^^ 
for low virtualities . The high values of the Cambridge algorithm, corresponds to large y -values, 
where the gluon jet is not the softest one. The CE and Boost algorithms both give results 
independent of y which are in very good agreement with O.BN^^ also for higher transverse 
momenta, b) Multiplicities in gluon jets as a function ofy for fixed k. The Durham algorithm 
assigns coherently emitted particles to the gluon jet, and the multiplicity increases with y. 
The Cambridge result is more independent of y, as long as the gluon jet is the softest one. 
of the Boost algorithms is very close to 0.5N'j^ for all y. 



All generated events are considered in the MC analysis. Neutrinos are excluded from the 
event, all other neutral particles are treated as massless while all charged particles are treated 
as pions. The obtained visible system in each event is boosted to its CMS before the analysis. 

To tag the gluon jet, we note that the algorithms are applicable to a general jet topology. It 
is therefore possible to study events where one jet is much softer than the other two, when 
it is a good approximation to assume the softest jet to be the gluon jet. This makes the 
analysis simple and independent of sophisticated tagging methods. The gluon jet is softest 
in the phase space region /t + 2||/| < L — ln(4). In our "soft tag" analysis, wc have restricted 
the phase space to k, + 2\y\ < L — 21n(4), in order to avoid events with two similar soft jets. 
The restriction still allows us to study scales up to k± ~ 20GeV. 

In an experimental situation, harder gluon jets can be identified using heavy quark infor- 
mation, why it is of interest to test the performance of the presented algorithms in a larger 
part of phase space. We have therefore also performed an analysis where the gluon jet is 
tagged using information available in the MC simulation, which however is experimentally 
non-observable. In this "angle tag" procedure, the jets are identified with the partons in 
such a way that the sum of jet-parton angles are minimized. 
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Figure 12: Comparison of jet identification 
methods. The cluster algorithms presented 
are designed to properly construct gluon 
jets in a large part of phase space. As- 
suming the gluon jet to be softest (symbols) 
restricts the available phase space to some 
extent, but the result is in very good agree- 
ment with more sophisticated gluon tagging 
methods (shaded areas). 



Figure 13: Multiplicities in the hardest 
quark jet obtained by the Cone Exclusion 
method compared to half an e~^e~ annihila- 
tion event. Results are plotted for several 
values of y. The range of results is given 
by a shaded region for the "angle tag". A 
reliable analysis can be performed in the 
kinematical region where the gluon jet is 
predominantly the softest one. 



Multiplicities in Gluon Jets 



In Fig |TT]a, the multiplicity of the gluon jet obtained by different algorithms is presented 
as a function of n. For each k, results for several values of y are plotted. The gluon jet 
was identified with the "angle tag" method. The Angular Ordering algorithms (Cambridge 
and Mercedes) have similar behaviours. They perform well at low transverse momenta, but 
start to show a y-dependence at higher virtualities. This is especially so for the Mercedes 
algorithm. The OSC algorithms (CE and Boost) give results independent of y and in very 
good agreement with O.bN^^. 



In Fig |Tl|b, the multiplicity of the gluon jet for fixed k is presented as a function of gluon jet 
rapidity y. To allow the analysis to include large y, the "angle tag" method is used. With 
the Boost algorithm, iVg (k) is independent of y and very close to the predicted |iVgg(K). 
The result from the Cambridge algorithm is independent of ?/ in a large range, but somewhat 
larger than ^N^. The steep rise of the Cambridge multiplicity at ?/ ~ 1.6 reflects the different 
treatment of the gluon jet when it is not the softest one. A conventional cluster algorithm - 
in this case the Durham algorithm - assigns particles from the region of coherent emission to 
the gluon jet. For large y this region becomes larger and the Durham multiplicity increases. 



In Fig O it is shown how the results using the very simple "soft tag" identification of gluon 



jets is in very good agreement with more sophisticated methods. 
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Figure 14: MC results of the multiplicity dependence on the jet resolution scale Kr for e~^e~ 
annihilation events with fixed CMS energy y/s = Mz- The gluon jet is identified both using 
the direction of the hardest gluon in the cascade and simply assuming the gluon jet to be 
softer than the quark jets. Solid lines are obtained using Eqs and (^^. a) The mean 
multiplicity in all events with no jet above Kr b) The mean multiplicity in all events with the 
hardest gluon jet at Kr- 



Multiplicities in Quark Jets 



In the Cone Exclusion algorithm, the OSC is used to define one-scale dependent regions 
also for the quark and antiquark jet. The results for the hardest quark jet are presented in 
Fig|l^. Again we note that jet mis- identifications in the "soft tag" method have little effect, 
and that this simple tagging procedure gives reliable results in a large kinematical region. 



7.3 Two-Scale Dependence 



The analysis required to compare data from fixed energy experiments with model predic- 
tions of the virtuality scale dependence is simple and straightforward. In Fig O the mean 



multiplicity in events where no gluon jet is found using a resolution scale /tr, and in events 
where the hardest gluon jet is found at Kr are presented. Since the full multiplicity of a 
two-jet or three-jet event is independent of how the particles are distributed among the jets, 
any fc^-based algorithm may be used. There is a fair agreement between MC simulations 
using the Durham algorithm, and the expectations from Eqs ( p5D and (26). 



With the Mercedes and Boost algorithms, the multiplicities in quark jets are expected to 
be 0.5iVqq(L -|- 2?/, k < Kr), where y and Kr are the kinematical variables of the gluon jet. 
Thus the energy evolution of the two-scale dependence can be studied in a fixed-energy 
experiment. The expected linear dependence in y (c.f. Eq (^)) is seen in Fig |15|, especially 
for the OSC-based Boost algorithm. 
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Figure 15: With the Mercedes and Boost algorithms, the energy scales of the quark jets 
are well defined and equals ln(s/A^) ± 2y, where the rapidity y of the gluon jet is defined by 
y = |ln[(l — Xq)/(1 — Xq)]. Thus it is possible to study the energy dependence of the two-scale 
dependent multiplicity. MLLA predicts this dependence to be linear, which also is the result 
for the Boost algorithm, while the Mercedes algorithm (based on AO) deviates somewhat at 
large \y\. 



8 Summary 



The ratio of the hadron multiphcity in quark and gluon jets is predicted to be 4/9 at very high 
energies. Corrections to this values calculated in the Modified Leading Log Approximation 
(MLLA) cannot expalin the difference between experimental data and this asymptotic value. 
In this paper we estimate the contribution from recoil effects to the ratio N^/Ng. Although 
formally of order 1/lns this contribution is quantitatively sizeable. Combined with the MLLA 
expression the result agrees well with MC simulations and with preliminary experimental 
data for jet energies equal to Mz/2 [|T^, We also discuss gluon jet definitions by which the 
scale evolution of gluon jets can be studied, using data from a fixed energy e~^e~ experiment. 
These jet definitions thus should enable an extended analysis which would be interesting to 



combine with the results of |Tl|] 



Another important point in this paper is that multiplicities in jets depend on two scales, the 
energy and the virtuality or maximal allowed transverse momentum. We derive expressions 
for the two-scale dependence and note that it can be examined with simple methods: Using a 
fc_L-based cluster algorithm to construct exactly three jets, the multiplicities can be examined 
ClS db function of jet transverse momentum. 

We discuss the "One-Scale Criterion" (OSC), which states that a one-scale dependent jet will, 
in some Lorentz frame, correspond to one hemisphere of a two-parton event where the energy 
and transverse momentum scales coincide. The relevant scale for the jet is the energy of the 
corresponding hemisphere. Most of the commonly used cluster schemes of today concentrate 
on reproducing the jet energy and direction, and put less emphasis on the assignment of the 
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soft particles to different jets. One exception is the Cambridge algorithm [Q, which based on 
Angular Ordering (AO) arguments constructs jets with one-scale dependent multiphcities. 
We present a set of algorithms designed to construct one-scale dependent jets. These are 
based either on AO or explicitly on the OSC. We examine the algorithms by analysing 
MC-generated events and comparing the obtained multiplicities in the jets with complete 
MC-simulated events at corresponding energies. 

Our study shows that all the presented algorithms perform well, but that the OSC methods 
are somewhat better than the AO ones for quantitative analyses of multiplicities in jets. 
With the OSC algorithms the gluon jet properties depend on only one scale. The treatment 
of quark jets differ, however. The "Boost algorithm" is particularly suited for studies of the 
two-scale dependence of quark jets, while the "Cone Exclusion" algorithm is designed for a 
study of one-scale jets, where the virtuality coincides with the energy of the jet. 

We note that the algorithms can be used in a large kinematical region. We do not require 
the events to be Mercedes-like or Y-shaped. In events where one jet is significantly softer 
than the others, this jet is predominantly the gluon jet. Performing the analysis on these 
events makes it less important to identify the gluon jet via quark taggings. 
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